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We considered a mesoscopic superconductor/normal metal (S/N) structure in which the N reser- 
| voirs are maintained at different temperatures. It is shown that in the absence of current between 

the N reservoirs a voltage difference Vr arises between the superconducting and normal conductors. 
' The voltage Vt oscillates with increasing phase difference <p between the superconductors, and its 

magnitude does not depend on the small parameter (T/cf). 

It is well known that if the terminals of a normal conductor are maintained at different temperatures, then in the 
absence of a current a thermoelectric voltage (V em /) appears between the terminals. The magnitude of V em f is equal 
O ' to ci(T/ep)ST/e, where C\ is a factor of the order 1, tp is the Fermi energy and ST is the temperature difference (see 
*_h ■ for example 

In this paper we analyse the thermoelectric effect in the mesoscopic structure shown in Fig.l. We show that when 
a temperature difference exists between the normal (N) reservoirs, a voltage between normal and superconducting 

! t j circuits Vr appears. Unlike in the normal case, the magnitude of this voltage does not depend on a small parameter 
(T/cf), and also this voltage oscillates as the phase difference (p between the superconductors varies. We assume that 
the superconductors are connected via a superconducting loop and the phase difference between them ip is controlled 
by an applied magnetic field. There is no current between the N reservoirs and the temperatures of the reservoirs are 
different: T(±L) — T Q ± ST. We will calculate the electric potential in the N film and, in particular, the potential Vr 

q . in the N reservoirs. Since we set the potential in the superconductors equal to zero, the potential Vr is the voltage 
difference between the N reservoirs and superconductors which arises in the presence of the temperature difference 
1 — '! ST. 

I In order to find the potential Vr, we need to determine the distribution functions f± and the condensate wave 
■ functions F R ^ induced in the N film. The distribution functions f± are related to the ordinary distribution functions 
for electrons n-\ and holes p[\ /+f = f+i = 1 — (n-\-\- pi) and /_f = = — (n-f — pi) (we assume that there is no 
spin-dependent interaction in the system). The function / + determines the condensate current and the function /_ 
determines the quasiparticle current and electrical potential (see for example Ref. 0] , where /+ and /_ are denoted by 
/ and fi respectively, and Ref. || where the application of the Green's function technique to the study of transport 
in S/N mesoscopic structures is discussed). These functions satisfy the kinetic equation (see Ref. [Q) 
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Ld x [M±d x f±(x)+JsMx) ± J an dMx)]= r[ A± 8{x - L x )+ A± 8{x + LJ]. (1) 



where all the coefficients are expressed in terms of the retarded (advanced) Green's functions: G R ( A ^ = G RiyA ^a z + 
pR{A). M± = {l _ qRqA t (j?flj?A )i)/2; j m = { F R F A ) z /2, J s = {\/2){F R d x F R - F A X F A ) Z , A ± - (vv s + 

9i?)(f±-fs±)-(9*kfs?+9*?U); 9i± = (l/4)[(F«±^)(F#±F#)] i; g z± = (l/4)[(F fl T F A )(F* ± F R )] Z . The 
parameter r = R/Rb is the ratio of the resistance of the N wire R and S/N interface resistance R^] the functions 
• • ■ A- and A+ coincide with A-,A+ if we make a substitution ip —> —p. We introduced above the following notations 
# > ■ (F R F A )i = Tr{F R F A ), {F R F A ) Z = Tr{a z F R F A ) etc.; u, v s are the density-of states in the N film at x = L x 
and in the superconductors. The functions fs± are the distribution functions in the superconductors which are 
assumed to have equilibrium forms. This means that fs+ = f e q = tanh(e/3 Q ) and fs- = 0, because we set the 
potential of the superconductors equal to zero (no branch imbalance in the superconductors). At the reservoirs the 
distribution functions f± obey the boundary conditions: f±(±L) = Fv±(±L), where Fy±(L) — (l/2)[tanh(/3 (e + 
eVr)) ± tanh(/3 (e - eV T ))],0 o = (2T o y\ 

In order to clarify the physical meaning of different terms in Eq.(l), consider the equation for /_. The first term in 
this equation is the partial quasiparticle current (the quasiparticle current at a given energy) in the N film . The second 
term is the condensate current and the third term is also the condensate current which appears under noncquilibrium 
conditions. The factors A- and A- are the partial currents through the S/N interfaces. The first term in A- is the 
partial quasiparticle current above (the term vvsf-) and below (the term gv^f-) the gap A, and the last term is the 
condensate current. The factors A + and A + are equal to zero below the gap (complete Andreev reflection). 
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In this paper we consider the case when the S/N interface resistance is greater than or approximately equal to the 
resistance of the N wire (r < 1). However r should not be too small because in Eq.(l) we neglected the inelastic 
collision integral. This implies that the condition 

{er,)- 1 « r < 1 (2) 

must be fulfilled; here r e is the energy relaxation time, e « min{T, and = D/L 2 is the Thouless energy. 
We consider the most interesting case of low temperatures (T « A) when for characteristic energies e << A the 
functions A + and A + equal zero as in this case vs ~ and F R = F A . The first integration of Eq.(l) yields 

M+d x f+(x)+J s f-(x) + J an d x f-{x)= J+ (3) 

where J+ is a constant of integration. In the limit of the small parameter r, all the corrections related to the proximity 
effect are small (the functions are proportional to r). Therefore in the main approximation in the parameter 

r, we have for the distribution function f+(x) 

f+(x) = 60e(x/L) cosh- 2 (/? e) + f eq (4) 

In obtaining Eq. (4) we have taken into account that in the N reservoirs the function / + (±L) has an equilibrium form 
with different temperatures:/+(±L) = f eq ± S(3e cosfT 2 (/3 e), where 6/3/0 = -6T/T,0 O = {2T )- X . Eq.(4) implies 
that the temperature gradient leads to a flow of nonequilibrium electrons and holes in the N wire. In order to find 
the distribution function f-(x), we integrate Eq.(l) 



M_d x f_(x) +J s f+(x)-J an d x f + (x)= Ji6(Li- \x\) + J 2 +,2-Q(\ x | -Li) (5) 

The constants J\ and J2+,2- are the partial currents in regions (— Li,+L\) , (L\,L) and (— L, — L{) respectively. 
Outside the interval {—L\, +L\) the " supercurrent" Js is zero (this follows directly from the expressions for J$ and 
for F R ). As there is no current between the N reservoirs, the integrals over energy e from J2+,2- should be equal to 
zero. In the absence of the temperature gradient we obtain from Eq.(5) /_ = and Jsfeq = Jieq = — r (9z- +9z+)feq- 
According Eq.(l) the constants J\ and J2+,2- are related to each other ( Kirchoff's law) 

J2+.2- - SJr = ±r[g 1+ f- - g z +5.f+](± Ll) (6) 

where 8J\ = J\ — Js and Sf + = f + — f eq . Integrating Eq.(5) outside the interval (— L\, +L\) and taking into account 
the boundary condition for f-(x) , we obtain J2+,2- = ±[/-(±L) — /-(±Li)]/i 2 - Using this expression and Eq.(6), 
one can show that in the main approximation in r the function /_ (x) is almost constant along the N wire and equal 
to /_(x) « /_(0) w f-(±L) = eV T p cosh" 2 (0 o e). Integrating Eq.(6) over energies we readily find eVr 

L ST roo 
e V T = -±— de <?2 + ecosfr 2 (/3 e)/ / de g 1+ cosh _2 (/3 e) (7) 
L ^ Jo Jo 

The potential Vr equals approximately the voltage difference between the N reservoirs and superconducting loop. 
It is worth noting that Vt, determined by Eq.(7), does not depend on the small parameter r because both functions 
g z+ and g 1+ are proportinal to r (however one should have in mind that according to the condition (2) this parameter 
must not be too small). The integrand in Eq.(7) can be calculated if the function F R is known from an approximate 
or numerical solution of the Usadel equation. In the limit considered, of small r, the retarded (advanced) Green's 
functions are readily found from the linearized Usadel equation. In this case we find 

g z+ =rRe(F* ^^ ) sin W g 1+ = r Im{F«[sinh(0 2 + 20,) + cosy sin hg 2 ] ^^ } (8) 

where 9 = 6\ + i#2,#i(2) = ^e^i(2) ; k e = y/— 2ie/eL,L 2 = L — L\ ,F R = A/ y/{e + iT) 2 — A 2 is the retarded Green's 
function in the superconductor. One can see that the ratio in Eq.(7) indeed does not depend on r. Numerical analysis 
of the Usadel equation shows that at a characteristic energy e = €l the difference between the linearized and exact 
numerical solutions of the Usadel equation is less than 10% even for r w 1. 

It follows from Eqs.(7-8) that the voltage Vt caused by the temperature gradient is zero when the phase difference 
between the superconductors is zero and oscillates with increasing ip. One can easily estimates the order of magnitude 
of V T . We find 
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eV T = STiLJL) sirup 



{ 



(T/e L )C 1 (p),T«e L 
(e L /T) 3 / 2 C 2 (ip),T»e L 



(9) 



Here Ci^if) are periodic functions of the phase difference ip of order 1 and are not zero when (p = 0. 
If we define the thermoemf Vt as the voltage between the N reservoirs and the superconducting circuit, we can state 
that this thermoemf is much larger than the thermoemf between the N reservoirs in the absence of superconductors 
because Vr does not contain the small parameter (T/ep) as is the case in a normal system (see for example Ref. 
HI). In addition, Vt oscillates with an applied magnetic field H {ip is proportional to H) allowing one to detect 
small temperature gradients. We stress once again that the thermoemf analysed in this work arises not between 
the normal reservoirs, as takes place in case of the ordinary thermoelectric effect, but between the superconducting 
and normal circuits. In Fig. 2 we plot the dependence Vjv on ip for various = (2T). In Fig. 3 the temperature 
dependence Vt at ip = n/2 is presented. We see that this dependence is non-monotonic with a maximum at T w cl 
(reentrant behaviour). We note that the influence of the proximity effect on the ordinary thermoelectric effect was 
studied theoretically in || . We ignore this effect regarding the thermoelectric current aVT as negligable, compared 
to the effect under consideration. In Ref. || the thermoelectric voltage was measured in complicated S/N structures 
which differ from the simple structure considered by us. The reentrant behaviour of Vt and its oscillations on ip were 
observed in this work. It is possible that the observed effects are related to those considered in this paper. It is 
worthwhile noting that the influence of the ordinary thermoelectric current aVT on the Josephson effect was studied 



The physical explanation of the effect is the following. The temperature gradient creates a deviation of the distrub- 
tion function 5f + = —(Sn + Sp) from the equilibrium form. On the other hand the superconductors do not affect this 
function because complete Andreev reflection conserves the total number of excess electrons and holes. The function 
<5/+ changes the condensate current flowing across the S/N interface. If <5/+ is created by an electrical current flowing 
between the N reserviors, it has the same sign at x = ±L and leads to a change in the Josephson current. In the case 
considered here the function Sf + has different signs at these points and leads to a variation of the condesate current 
of the same sign at different S/N interfaces. Therefore, the potential Vt arises in the N wire producing a subgap 
current rt/i + /_(±Li) which compensates the current SJs- 

In summary, we have calculated the voltage Vt between the superconducting and normal curcuits in a S/N meso- 
scopic structure where the normal reservoirs are maintained at different temperatures. This voltage arises due to a 
branch imbalance in the N film and oscillates with varying phase difference. Its magnitude does not contain the small 
parameter (T/ejr) which is present in normal systems and is of the order 5T(L\/ L) / e. 

We are grateful to the Royal Society and to the EPSRC for their financial support. We are grateful to Yu.Gal'perin 
for his useful comments and to V.V.Pavlovskii for his assistance. 



[1] A.A.Abrikosov, Fundamentals of the theory of metals, North-Holland, 1988 

[2] A.I. Larkin and Yu.N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 47, 1136 (1964) [Sov. Phys. JETP 20, 762 (1965) 

[3] C.J.Lambert and R.Raimondi, J. Phys., Condens. Matter, 10, 901 (1998) 

[4] R.Seviour and A.F.Volkov, Phys.Rev 61, R9273 (2000) 

[5] N.R.Claughton and C.J.Lambert, Phys. Rev. B53, 6605 (1996). 

[6] J. Eom, Chen- Jung Chien, and V.Chandrasekhar, Phys. Rev. Lett. 81, 437 (1998) 

[7] A.G.Aronov and Yu.M.Gal'perin, JETP Letters 19, 165 (1974); M.V.Karstovnik, V.V.Ryazanov and V.V.Shmidt, JETP 
Leteers 33, 357 (1981) 



long ago @. 



3 



V T s 



s v T 



T0-8T 



To + 5T 



FIG. [i] Schematic view of the 4-terminal S/N/S structure under consideration. The electric potential of the superconductors 
is zero. The N reservoirs are disconnected from the external circuit. 




FIG. ^The dependence of the normalised thermoeletrical voltage Vt = eVr / {ST L\ / L) on the phase difference ip for various 
(3 — T/tL (the parameters are A/el = 10, Li/L — 0.5, r = 0.3). 




FIG. ^| The temperature dependence of the normalised voltage Vt at ip = ir/2, for the same parameters as in Fig. 2. 
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